Phase diagram of the Holstein polaron in one dimension 
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The behavior of the ID Holstein polaron is described, with emphasis on lattice coarsening effects, 
by distinguishing between adiabatic and nonadiabatic contributions to the local correlations and dis- 
persion properties. The original and unifying systematization of the crossovers between the different 
polaron behaviors, usually considered in the literature, is obtained in terms of quantum to classical, 
weak coupling to strong coupling, adiabatic to nonadiabatic, itinerant to self-trapped polarons and 
large to small polarons. It is argued that the relationship between various aspects of polaron states 
can be specified by five regimes: the weak-coupling regime, the regime of large adiabatic polarons, 
the regime of small adiabatic polarons, the regime of small nonadiabatic (Lang-Firsov) polarons, 
and the transitory regime of small pinned polarons for which the adiabatic and nonadiabatic contri- 
butions are inextricably mixed in the polaron dispersion properties. The crossovers between these 
five regimes are positioned in the parameter space of the Holstein Hamiltonian. 



I. INTRODUCTION 

Recent years have witnessed a constant growth of at- 
tention for polaron- related physics. These investigations 
are motivated by many examples of experimental evi- 
dence of polaron formation, found for broad classes of 
compounds like quasi-lD conductors (MX-chains, con- 
ducting polymers), ionic crystals, transition-metal ox- 
ides, and fullerenes. Furthermore, significant isotope ef- 
fects in the behavior of charge carriers are observed for 
magnetoresistive manganites 1 -^ and high-temperature 
superconducting cuprates 4 ^^, which indicates the pres- 
ence of polaronic effects involving spatially localized cor- 
relations between the charge carriers and the lattice 
deformation field. Although the polaronic correlations 
alone are usually not believed to be sufficient to ac- 
count for colossal magnetoresistance 7 - or for the high-T c 
superconductivity^^, their investigations are neverthe- 
less of particular interest in this context. 

Depending on the model parameters, the polaron con- 
cept itself involves very different types of behavior, which 
complicates the analysis of experimental results. As one 
would expect, the polaronic features are the most promi- 
nent for strong electron-phonon couplings. In particular, 
in this limit, the polaron theory predicts a broad contri- 
bution to the mid-infrared optical conductivity, accom- 
panied by a strong suppression of the Drude pea k 10 ' 11 . 
Such behavior is indeed observed for various compounds. 
For example, the formation of small and large polarons 
has been reported 1 ^ for Lai-^Ca^MnOa (LCMO) and 
Lai-^Sr^MnOs (LSMO), respectively. Large polarons 
have been invoked for Nb doped SrTiOy^. The simul- 
taneous presence of large and small polarons has been 
argued for some oxides like Pr 2 Ni04.2z^- The polaronic 
nature of the charge carriers has also been reported for 
slight doping from optical measurements of the high-T c 
cuprate Nd2-xCe x Cu04_ J/ (NCCO) 1 ^. Strong ionic cou- 
plings do occur in ionocovalent materials such as high-T c 
superconductors 1 ^ 

In the context of ARPES measurements, small po- 



larons have been invoked to account for the quasi-particle 
peak with very narrow dispersion near the Fermi level 
and the broad peak that roughly follows the free-electron 
ban d 17 ' 18 , observed in the undoped cuprate*^. Small po- 
laron behavior involving very strong couplings has been 
found for the undoped Ca2Cu02Cl2 cuprate^ -, the quasi- 
one-dimensional conductor (TaSe^z^ and molybdenum 
bronze^ 2 -. In the context of neutron and X-ray measure- 
ments, the observed diffuse scattering can be related to 
the polarons involving local lattice distortions that are 
slow on the time scale of typical phonon vibrations, as 
suggested for colossal magnetoresistive manganite o 23 ' 24 . 
In the case of conducting polymers, the lattice defor- 
mation and the softening of local phonon modes as- 
sociated with the solitons and, possibly, polarons has 
been deduced from the doping-induced infrared-active 
lattice modes in photoinduced absorption spectr a 25 ' 26 . 
Recently^ 7 ., using the femtosecond optical spectroscopy 
the internal dynamics of the polarons in GaAs, charac- 
terized by weak electron-phonon interaction, was inves- 
tigated. This technique seems promising for studies of 
polarons in other compounds too. 

On the theoretical side, many ongoing inves- 
tigations are focused on the Holstein polaron 
problem 2 ^^^^^ 3 ^. 3 ^ 7 .. Although described by 
only two independent parameters, which necessarily sim- 
plifies the behavior in comparison with real systems, the 
Holstein model contains the essential physics governing 
the polaron formation in the presence of short-ranged 
ionic coupling between the electron and optical phonons. 
In particular, it uses the band picture for the electron, 
while the discrete lattice deformation field is treated 
quantum mechanically. As a result, one finds that the 
Holstein polaron exhibits fundamentally different behav- 
iors, ranging from quantum to classica l 38 ' 39 , from weak 
coupling to strong couplin g 40 ' 41 ' 42 ' 43 ' 44 , from adiabatic 
to nonadiabati c - 40 ' 41 ' 42 ' 43 , from itinerant to self-trapped 
polaron a 42 ' 45 , and from large to small polaron a 42 ' 45 ' 46 . 
The present work provides a unifying description of all 
these aspects of polaron physics within the corresponding 
phase diagram. In order to improve the understanding of 
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the interrelation among various parameter regimes and 
to situate the crossovers between them, many known 
results are reviewed and supplemented with additional 
details. In this way the applicability of various approxi- 
mations to particular regimes is fully clarified. The key 
elements of the formation of the Holstein polaron are 
identified, resulting in an original and comprehensive 
interpretation of the low-frequency polaron dynamics. 
The features of the polaron band structure are explained 
so in detail. 

Two perturbative approaches are commonly applied in 
the context of the Holstein polaron problem. The first 
involves the expansion around the atomic limit, while 
the second treats the electron-phonon coupling g as a 
perturbation in the theory, which is translationally in- 
variant from the outset. Although this latter theory de- 
scribes the nonadiabatic-adiabatic crossover in the con- 
tinuum limit (large Holstein polaron) appropriately and 
indicates the continuous to discrete adiabatic crossover 
correctly, the actual calculation of the lattice coarsening 
effects proved too intricate^. On the other hand, when 
starting from the atomic limit it is difficult to reach the 
adiabatic regimes, in particular, it is difficult to describe 
adiabatic contributions to the dispersion even for small 
polarons. 

Various approximate diagrammatic techniques have 
been developed to account for high order corrections in g, 
including the self-consistent Born approximation 4 ^, the 
dynamical mean field theory 4 ^, and the momentum aver- 
aging approac h 48 ' 50 . These approaches have in common 
that the electron self-energy is treated as a local quantity. 
However, as shown in Refs i 47 i 51 , it is necessary to retain 
the non-local contributions to the electron self-energy in 
order to correctly describe the adiabatic polaronic corre- 
lations involving multiple lattice sites. In fact, with the 
increasing range of adiabatic correlations, many numer- 
ical methods proposed in the literature become inade- 
quate, including those based on selective diagrammatic 
calculations, as well as those based on direct calculations 
of polaron states. The well controlled general way of 
dealing with this problem was recently proposed in the 
context of the relevant coherent state method (RCSM) 5 -^, 
which calculates accurately the low-frequency polaron 
spectra and wave functions for the whole parameter space 
of the Holstein polaron. For this reason, the RCSM is 
used in the present study as the numerical tool for fulfill- 
ing the gaps in the phase diagram that are not covered 
by the limiting analytical solutions. 

Until now, the various regimes of the Holstein polaron 
have not been described on an equal footing. By review- 
ing the results for some of them and developing further 
the analysis for others, we are able to achieve a compara- 
ble level of understanding for all regimes across the entire 
parameter space. Some entirely new results for the small 
adiabatic polaron and the crossover between its adiabatic 
and nonadiabatic translational dynamics are presented 
here. Regarding large adiabatic polarons, the formalism 
of the moving set of coordinates is used to derive the 



expression for the polaron band structure in presence of 
lattice coarsening effects, an important problem also ap- 
pearing in the soliton theory. Furthermore, the nature of 
singularities appearing within the formalism of the mov- 
ing set of coordinates in the large and in the small polaron 
regime is clearly distinguished. 

In Sec. Hi] the model and some general properties of 
the Holstein polaron states are discussed, whereas a 
short overview of the perturbative results is presented in 
Sec. IIIII The key-elements that characterize the adiabatic 
polaron dynamics are identified in Sec. IIVI The large 
adiabatic polaron and the two crossovers, towards the 
weak-coupling regime and the regime of small adiabatic 
polarons, are analyzed in Sec. E] In Sec. IVIl the small 
pinned polarons are investigated in detail, with particular 
emphasis on the characteristics of the polaron dispersion 
properties. The RCSM results are used to disentangle 
various contributions to the polaron dynamics along the 
crossover between freely moving and nearly pinned po- 
laron states, which spans across the phase diagram from 
the nonadiabatic to the adiabatic limit. The overview 
of the most important results and conclusions is given 
with the phase diagram for the one-dimensional Holstein 
polaron problem in Sec. IVII} where the generalization 
of the phase diagram to higher dimensional cases is also 
discussed. A short summary is presented in Sec. IVIIII 

II. GENERAL 

The present work considers polaron properties within 
the Holstein model^ in one dimension. In momentum 
space, the Holstein Hamiltonian is given by 

H = -2t^cos(fc) c\c k + Lu ^b\b q 

k q 

- s/Viv^cl+^C^ + M. (i) 

where N is the number of lattice sites (N — > oo). c\ 
and are the creation operator for the electron and 
phonon, respectively. Only the single-electron problem 
is considered. The electron is described in the tight- 
binding nearest-neighbor approximation, with t the elec- 
tron hopping energy. The coupling with the branch of 
dispersionless optical phonons is local, g is the electron- 
phonon coupling energy. The phonon energy ujo is taken 
throughout this work as the unit of energy. Besides t/uo 
and g/ujQ, the two other parameters used are the bind- 
ing energy of the polaron in the atomic limit (t = 0), 
e P = 5 2 / w o, and the dimensioneless parameter A = e p /t. 
Unlike g, e v and A are adiabatic parameters in the sense 
that they do not depend of the mass of the nucleus. 

The part of the spectrum below the phonon threshold 5 -^ 
for incoherent scattering, ui < E + lu , involves eigen- 
states for which the phonons are coherently correlated 
with the electron in the space. In other words, the 
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electron-phonon correlation length d is for these eigen- 
states finite. This lenght defines the spreading of the 
(polaronic) phonon field that moves with the electron. 
In particular, the spectrum for ui < Eq + luq can be inter- 
preted in terms of one (or more) polaron bands^. Each 
polaron state in the band is an eigenstate of the total mo- 
mentum K = J2k k cl c k + J2 q 9 b\b q , where K commutes 
with the Hamiltonian {TJ. 

Since the homogenous q = phonon mode in Eq. |TJ, 
b q=Q = Y^ n ^n/y/N, interacts only with the total electron 

density that is conserved, X)fe c l c fc = 1> the correspond- 
ing part of the Hamiltonian 

H g=0 = lu bl= b g=0 + g [b\ =0 + b q=0 )/VN , (2) 

can be solved exactly^. The solutions are given by the 
ground and excited states of the displaced q = har- 
monic oscillator. The displacement of the q = oscillator 
defines the mean total lattice deformation, 

xtot = x Vn (*|(^ =0 + b g=0 )\$>) = 2.t g/ujQ , (3) 

where \^} is any eigenstate of the total Hamiltonian (Q]) 
a nd xq is t he root mean-square zero-point displacement 
-\/l/2Af(jJo, with M the mass of the nuclei (h = 1). 

In the present work, the local interplay between the 
electron and its accompanying phonon field is referred 
to as the local dynamics. In the static u> — > Eq limit, 
the local electron-phonon correlations are characterized 
by the mean total lattice deformation x to t and the cor- 
relation length d. While xtot is simply determined by 
the ratio g/ujQ, the correlation length behaves differently 
in the case of the adiabatic and the nonadiabatic local 
correlations. The polaron dynamics is adiabatic when 
the electron follows the motion of the lattice deformation 
almost instantaneously. Nonadiabatic dynamics involve 
processes for which the electron is temporarily detached 
from the deformation. 

The distinction between adiabatic and nonadiabatic 
contributions is important for the understanding of po- 
laron dispersion properties as well. By analogy to d, 
which characterizes the local correlations, the effective 
polaron mass m po i describes the translational dynamics 
involving the joint motion of the electron and the phonon 
field along the lattice. Generally, the effective mass is 
given by a different combination of parameters in the adi- 
abatic and nonadiabatic cases. In addition, the polaron 
dispersion properties depend on the value of the electron- 
phonon correlation length d (polaron size). Namely, po- 
larons that are large with respect to the lattice constant 
a move as free particles and can be investigated in the 
continuum approximation. However, lattice coarsening 
effects come into play for small polarons d ~ a, because 
they can be pinned to the lattice. 

The polaron formation is discussed in this work in 
terms of analytical results for the ground state energy, the 



polaron size, and the effective mass for various parame- 
ter regimes. However, in the crossover regimes, rather 
than simply focusing on these quantities, the analysis 
is performed by considering the polaron band structure 
for lj < Eq + u>o. In this way, important insights into 
the nontrivial mixing of fundamentally different contri- 
butions to the polaron dynamics, which would not be 
apparent from an analysis of the ground state and the 
effective mass alone, are obtained. 

III. PERTURB ATIVE APPROACHES 

In order to discuss the nature of polaron states in the 
weak-coupling and small t limits, it is convenient to ana- 
lyze the leading corrections to the polaron binding energy 
and the effective mass in the context of two standard per- 
turbative approaches. In particular, the perturbative cor- 
rections in g are discussed in Sec. IIII A\ while the power 
series expansion in t is examined in Sec. IIITBl In the first 
case, the perturbative theory is translationally invariant 
from the outset, while in the second case the perturbative 
expansion starts with the atomic t = limit, i.e., from 
the state which breaks the translational symmetry. In 
this connection, it is interesting to point out that in the 
small g and t regime (g,t <C Wo) the leading corrections 
for both perturbative approaches give the same polaron 
behavior. This result serves as good evidence support- 
ing the view that both perturbative series are equivalent, 
when summed to infinite order. 



A. Leading corrections in g 

Up to the second order in g of the nondegenerate 
Schrodinger perturbation theory^, the polaron binding 
energy is given by££ 

A po i = e p /y/l + 4t/Lj , (4) 
while for the effective polaron mass one obtains 

nipoi 2t + uj Q 

—^— = l + e p j (5) 

m e i (cu )2 (At + u) )2 

with £ p = g 2 /u>o. Up to the second order, the electron- 
phonon correlation length d is independent of the cou- 
pling g. In particular, in the static u> — > Eq limit, the 
mean lattice deformation decays exponentially with the 
distance m from the electron^, 



u m ~ ^l + wo/2t- \JujQ/t + ul/±t 2 j . (6) 

The polaron is small (d/a 1) for t/ui <C 1, and large 
(dja > 1) for t/ujQ > 1. 
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In the small polaron limit t/u>o <C 1, the binding energy 
and the effective mass behave according to 



A ^ _ m pol ^ -, , A po ; 
<-*pol ~ C'p i ~ - 1 - T" 



Tig/ 



UJ 



(7) 



In the opposite, large polaron limit d/a ~ ^Jtjujo 3> 1, 
one obtains 



m po i 
m e i 



1 + 1 A p°' 
2 w 



(8) 



By introducing iV^ = 1 + 2<i/a, where Nd is the ef- 
fective number of lattice sites involved by the electron 
dressing effects (the polaron size) , the binding energy for 
weak couplings can be, to a good approximation, written 
as Ap i f=a e p /Nd, satisfying the both limiting behaviors 
given by Eqs. |f7| and ([8]). The fact that the correlation 
length d depends on the lattice mass through the ra- 
tio t/uiQ reflects the nonadiabatic nature of the electron- 
phonon correlations in the weak-coupling limit. This 
holds for any dimensionality of the system. In partic- 
ular, in the context of the diagrammatic perturbation 
theory, one finds that the leading contribution to the 
electron self-energy is given by the local bare electron 
propagator—, T,^(uj) = g 2 Gi oc (uj— uj a ). In the weak cou- 
pling limit, the ground-state properties are determined 
by the value of £( 2 )(w) at the bottom of the nearly free 
electron band, which leads to the correlation length d 
independent of g^l. 

For large polarons, Eq. |(8j) can be derived in the con- 
tinuum approximation by considering the long- wave limit 
for the electron dispersion^. On the other hand, for 
small polarons, the band structure of Eq. (jTJ should be 
preserved in order to derive Eq. ([7])- 

The weak-coupling regime, for which the electron is 
only weakly dressed by the phonon cloud, is determined 
by the condition A po i < ujq. For the one dimensional 
system, one obtains 



ep/ui Q < 7V d « 1 + 



UlQ 



(9) 



i.e., g/uj < 1 for t <C ui and g/uj < {t/^o) 1 ^ 4 for 
t^>uj . Namely, for A po ; ~ uj , the second and fourth or- 
der contributions to the binding energy are comparable^, 
meaning that the two-phonon processes are of similar im- 
portance as the single-phonon processes. In other words, 
for Ap i > uiq, the electron is strongly renormalized by 
the phonon field. 



Eq = —Bp. By taking into account the first order correc- 
tions in t, one obtains the well-known nonadiabatic Lang- 
Firsov small polaron band. With respect to the free- 
electron case, the bandwidth of the Lang-Firsov polaron 
band is reduced by the electron quasi-particle weight, 



t LF = t CXp (-e p /u ) 



(10) 



The polaron hopping in Eq. ifTOj) involves processes for 
which the electron hops to the neighboring site by de- 
taching itself from the phonon field. Such nonadiabatic 
contributions delocalize the whole polaron, rather than 
breaking the coherent correlation between the electron 
and phonons, since below the phonon threshold the inco- 
herent scattering is forbidden by the energy constraint. 

In the small t and weak coupling e p <C u>q limit, the 
first order expansion in t computed by the degenerate 
Schrodinger perturbation theory gives the same polaron 
binding energy and the effective mass (i.e., by expand- 
ing Eq. flOj) ) as in Eq. |f7|, derived in the context of 
the perturbation theory in g. The crossover between the 
weak-coupling regime and the strong-coupling regime, in 
which the (Lang-Firsov) polaron is characterized by an 
exponentially reduced bandwidth, occurs for 



9 ~ • (11) 

In the Sp/ujQ ^> 1 limit, the polaron hopping energy 
tLF, given by Eq. lfTD|). becomes negligible. That is, for 
£ p ^> u>q > t, the non-exponentially small correction to 
the polaron ground state energy is given by the second 
order contribution in t that is independent of the polaron 
momentum^, 



En 



t 2 

-e p -2^£[Ei (2e p /ujo) 
-e p (1 + t 2 /e 2 p ) . 



7o - In (2e p /w )] 



(12) 

with Ei(x) the exponential integral, and 70 Euler's con- 
stant. The second term in Eq. lfl2|) is independent of 
the lattice mass, revealing the adiabatic nature of local 
electron-phonon correlations. One notices that this be- 
havior is fundamentally different from that obtained in 
the weak coupling limit. Namely, in the latter case, t 
appears in the polaron binding energy in Eq. |(4]) through 
the nonadiabatic ratio t/ujo- 

The nonadiabatic hopping fi"0j) is apparently indepen- 
dent of the system dimensionality D. On the other hand, 
the second order corrections in Eq. lfl2|) are D-dependent. 
However, for e p S> t^o ^ t, their nature is adiabatic just 
like in Eq. (fl2| . irrespectively of the system dimension. 



B. Leading corrections in t 

In the atomic limit t = 0, the polaron formation 
involves a single lattice site and the ground state en- 
ergy of the electron-phonon system is TV-fold degenerate, 



IV. ADIABATIC ALLY SELF-TRAPPED 
ELECTRON 

When g and/or t cannot be treated as small perturba- 
tions, instead of dealing with the (high-order) perturba- 
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tion theory^, it is convenient to approach the polaron 
problem in terms of the adiabatic approximatio n 60 ! 61 . 
Namely, when the electron spectrum is calculated as 
a function of the lattice deformation, as illustrated in 
Fig. [H a large gap Aad opens for strong couplings. The 
gap Aad separates the excited states from the ground 
state, the latter being characterized by a finite electron- 
phonon correlation length d. This length is, unlike for 
weak couplings, independent of the lattice mass, d ~ 1/A, 
with A = E p /t = g 2 /t u> . 

In the adiabatic limit the electron always stays in the 
localized ground state, i.e., the electron is self-trapped. 
The nonadiabatic corrections involve transitions of the 
electron back and forth into the excited states of the adia- 
batic electron spectrum [see Fig.[T|. The adiabatic regime 
is frequently assigned to large values of the ratio t/ujQ 
(fast electron vs. slow lattice). However, it is important 
to realize that the condition Aad 3> wo that makes the 
nonadiabatic corrections small depends on the coupling 
g, with Aad ~ £p for small polarons and Aad ~ (£p A) 
for large polarons. This means that for small polarons 
Aad can be large even in the small t/ui limit. On the 
other hand, for large polarons, the condition t/u>o ^> 1 
by itself is not sufficient to guarantee Aad ^> uiq, i.e., 
the adiabatic behavior. 
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Figure 1: (Color online) The polaron formation in the Hol- 
stein model for strong-couplings. For all lattice deforma- 
tion relevant for the polaron dynamics the associated electron 
spectrum, calculated as a function of the lattice deformation, 
exhibits a large gap Aad ^> Wo, separating the ground (lo- 
calized) state from the excited states. Although the electron 
is trapped by the lattice deformation, the translational sym- 
metry is preserved because the electron can move along the 
lattice together with the lattice deformation. 



A. Adiabatic approximation 

Within the adiabatic approximation, the electron wave 
function is determined by the lattice deformation, while 
the dependence on the lattice momentum is assumed to 
be negligibly small. Accordingly, the wave function of 
the polaron is decomposed into the product of the lattice 
Q ph (u) and the electron rj n {u) part, 

=& h (u) ® 5>„(u)ct |0) , (13) 

n 

with r] n (u) the ground-state electron wave function. 
Hereafter, for sake of brevity the lattice deformation is 
denoted by a vector u = {«„}. That is, u represents 
a point in the A^-dimensional configuration space of lat- 
tice deformations. The component u n is the lattice de- 
formation at the site n expressed in dimensionless unit; 
u n = 1/2 corresponds to the root-mean square zero-point 
displacement of the free lattice. In present notation, |w| 
is the norm of the vector, |u| = (J2 n u n)^- The unit 
vector is denoted by u, \u\ = 1. 

The expectation value of the Hamiltonian ([1]) with re- 
spect to the electron ground state r\ n (shifted by the zero 
point energy) is given by 

n fl n 

- t^VniVn+i +Vn-i) - 2.9 ^ u n \rj n \ 2 (14) 

n n 



The requirement that the electron energy sad (u) is min- 
imal for the given lattice deformation u yields 

£ad{u) r] n = -t {r/n+i + r/n-i) - 2,g u n r\ n , (15) 

the solutions of which can easily be found numerically. 
The solution with the lowest energy defines the ground 
state, whereas other solutions define the excited states of 
the adiabatic electron spectrum, shown schematically in 
Fig.CQ 

Since the electron ground state r\ n in Eq. I|14p is a func- 
tion of the deformation u, with the exception of the first 
term, the rest of Eq. (JT4J) can be interpreted as the adia- 
batic potential Uad(u) that characterizes the lattice dy- 
namics in the adiabatic approximation. The Hamiltonian 
Had is translationally invariant just as is the Hamilto- 
nian |TJ. The first term in Eq. (fl4| represents the kinetic 
energy of the lattice. The contribution to Had arising 
form the fact that the electron wave function and the lat- 
tice kinetic energy operator do not commute is neglected 
in the adiabatic limit. 

For the Holstein model, the homogeneous (q = 0) 
lattice mode couples to the total electron density. For 
this reason one should consider just those lattice defor- 
mations in Eq. (fT5"ll that satisfy the sum rule ((3]), i.e., 
u n = g/wo- The wave functions r\ n depend only on 
A, while the electron spectrum exhibits a simple scaling 
property, 

SAD(u,e p ,X) = s p eAD(^ou/g, A) . (16) 



Furthermore, the same property is exhibited by the adi- 
abatic potential, 



Uad(u) = loq I u | 2 +e A d(u) = s p U A D(u u/g, A) . 

(17) 

The only parameter that governs Uad(u) is A, namely 
the parameter that defines the adiabatic electron-phonon 
correlation length (the polaron size). 



B. The moving set of coordinates 

In order to discuss the adiabatic polarons, it is con- 
venient to introduce a moving set of coordinates u = 
{s,Q Q (s)}, which allows a natural distinction between 
the translational and local dynamics. Namely, in this 
formalism, s represents the position of the polaron in 
its translational motion along the minimal energy path. 
Q a (s) are the coordinates of the normal modes that move 
with the polaron. These coordinates describe the local 
lattice dynamics involving displacements orthogonal to 
the minimal energy path associated with the translation. 

The formalism of the moving set of coordinates was 
originally introduced in the context of the quantum field 
theory in order to analyze soliton-like solutions involv- 
ing multi-dimensional configuration spaced. It is exten- 
sively used in quantum chemistry to calculate reaction 
rates for adiabatic potentials that describe various chem- 
ical reaction o 63 ' 64 . In these problems, the minimal energy 
path (reaction path) usually connects two minima, sepa- 
rated by a potential barrier, corresponding to the desired 
reactants and products. The particularity of the polaron 
problem is that the minimal energy path involves motion 
along a periodic lattice, 



Uad{s, Qa{s)) = U A d(s + S ai Q a (s + 8 a )) 



(18) 



with s denoting the length of the path in the con- 
figuration space u that corresponds to the translation 
of the polaron by one lattice site. The approach with 
the moving set of coordinates was used for the adi- 
abatic polaron problem previously in the continuum 



approximation 



65,66,07,08 



. The analysis presented here ex- 



tends those investigations by addressing the issue of lat- 
tice coarsening effects. 

Due to the periodicity, it is sufficient to analyze the 
properties of the adiabatic potential V ad for one unit cell 
in order to identify the key elements of the polaron adi- 
abatic dynamics. In particular, in the following sections 
three aspects of the adiabatic dynamics are separately 
discussed: the emergence of the barrier for the transla- 
tional motion due to the lattice coarsening, the phonon 
softening effects related to the local dynamics involving 
fluctuations around the minima of the adiabatic poten- 
tial, and the appearance of the kinematic coupling^ be- 
tween the translational and local dynamics due to the 
curvature of the minimal energy path. 



C. Peierls-Nabarro barrier 

The translational motion of the polaron along the min- 
imal energy path is characterized by iV-fold degenerate 
minima and saddle points of Uad(u). Their difference 
in energy defines the Peierls-Nabarro (PN) barrier Apjv, 
which is the minimal energy barrier that must be over- 
come classically in order to move the adiabatic polaron 
from one lattice site to another. By choosing one of the 
minima of Uad{u) as the origin for the coordinate s, the 
PN barrier in terms of the moving set of coordinates may 
be expressed as 



Apjv = U AD (s a /2, Q a = 0) - Uad(P, Q a = 0). (19) 

The lattice configurations associated with stationary 
points of the adiabatic potential are obtained by extrem- 
izing Uad{u) with respect to u n , 



VU AD {u) = Q 



U n = (ff/^o) \Vn 



(20) 



After inserting Eq. (|20p into Eq. lfl5j) one obtains the 
discrete nonlinear Schrodinger equation whose stationary 
solutions define the stationary points of Uad{u). The 
lattice deformations u M corresponding to the minima are 
centered at the lattice site, while the deformations u PN 
corresponding to the saddle points with the lowest energy 
are centered between the neighboring lattice sites. By 
inserting Eq. l(20|) into Eq. lfl5|) . one can easily verify 
that the PN barrier scales as 



A PN = U AD {u™)~ U A d(u m ) = £ P A PN {\) . (21) 

The PN barrier rescaled by e p in Fig. [2] is shown as a 
function of A. In the small polaron limit A> 1, the PN 




Figure 2: The PN barrier Apjv(A) = Apjv/e p as a function 
of A. It is exponentially small for A < 1, whereas for A 3> 1 it 
asymptotically approaches the constant value Apn ~ 1/2 
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barrier is the highest, being of the same order of magni- 
tude as the polaron binding energy, Ap« w e p /2. For 
large polarons, d/a ~ A -1 » 1, the shape of the lattice 
deformation is not significantly altered by the lattice dis- 
creteness as the polaron moves along the minimal energy 
path. Accordingly, in the small A limit one finds that the 
PN barrier is exponentially smal l 70 ' 71 , 

lim Apjv(A) - A- 2 cxp(-27r/A) . (22) 

It is worth noting that A < 1 appears also in the transla- 
tionally invariant diagrammatic theory as the condition 
for the validity of the continuum approximation in the 
adiabatic limit^ 7 -. 



D. Phonon softening 

Apart from the role in the polaron translational dy- 
namics, the adiabatic dynamical correlations manifest 
themselves through the softening of the local phonon 
modes associated with the lattice deformation. The adi- 
abatically self-trapped electron introduces an effective 
coupling between lattice vibrations within the range of 
the electron-phonon correlation length d ~ 1/A. As a 
consequence, the effective elastic constants which charac- 
terize the lattice-deformation normal modes are weaker 
than in the case of uncorrelated loq lattice vibrations. 

In the limit of the large PN barrier Apjy ^> Woj the 
local dynamics of the polaron is restricted to harmonic 
fluctuations around the minima of Uad{u), whereas an- 
harmonic contributions to the local dynamics responsible 
for the polaron derealization can be neglected. Namely, 
because of the exponentially small probability for the 
tunneling through the PN barrier, the time scale 1/wo 
associated with the local dynamics is much shorter than 
the time scale on which the polaron moves along the lat- 
tice. The harmonic expansion of the adiabatic potential 
Uad(u) around the minimum Uad(u m ) is given by 

Uad{u M + h) « U A d{u M ) +uj hKh. (23) 

h represents the displacement from the equilibrium con- 
figuration u M . K is the elastic constant matrix, 

K, j = 8 id + II,., . (24) 

Here, Sij is the Kronecker delta. Hij is the (static) elec- 
tron polarizability matrix describing the adiabatic cor- 
relations. It is determined by the adiabatic changes of 
the electron wave function in response to infinitely small 
lattice displacements from the equilibrium configuration 
u M . In terms of the electron basis corresponding to u M , 
one obtains^ 

n,,(A) = -4e p?7 ^£-^, (25) 




0.5 1 1.5 2 X 



Figure 3: (Color online) The eigenfrequencies ujp/uio of the 
six lowest normal modes (5 as functions of A. The pinning 
and the breather mode are denoted by P and B, respectively. 

where the index denotes the ground state solution of 
Eq. l(T5)l . whereas v > denotes the excited state solu- 
tions. Because the electron spectrum e , e v , scales ac- 
cording to Eq. l(T6)) . IT j (A) is a function of A. 

By solving the matrix problem ((24)) . one obtains the 
normal modes of the lattice deformation, 

U AD {Qp) « U A d{u m ) + Y,^Q}- (26) 



It should be stressed that the harmonic term in Eq. ((26)) 
includes the (local) fluctuations along the minimal energy 
path too. That is, the coordinate of the normal mode (3 
with the lowest frequency in Eq. ((26)1 corresponds to the 
coordinate s in the formalism of the moving set of co- 
ordinates {s,Q a (s)}. In this respect, it is important to 
distinguish between the fluctuations around the minima 
associated with the coordinate s and the motion from one 
minimum to another involving the same coordinate. The 
latter contributes to the translational dynamics, whereas 
the former is associated with the local dynamics. Here- 
after, the index /3 is used to denote the full set of N 
normal modes appearing in Eq. (f26)l . calculated with re- 
spect to the potential minimum. The index a involves 
the moving set of N — 1 normal modes orthogonal to the 
minimal energy path. 

The lattice deformation Qp = {Qp^n} associated to 
the normal coordinate Qp may be either symmetric or 
antisymmetric with respect to the polaron center, which 
defines the parity of the normal mode (3. Because the 
elastic constant matrix ((24)) is determined by A, the same 

applies to the unit vectors Qp and to the ratios uj^/uq 
defining the frequency softening. 

The eigenfrequencies ojp/wo of the six lowest normal 
modes are plotted in Fig. [31 The softening effects are 
particularly strong for the lowest odd (pinning, /3 = P) 
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and the lowest even (breather, (3 = B) mode. The fre- 
quency of the odd pinning mode u>p = p tends to zero in 
the large polaron limit [the left side of Fig. [3] . This is a 
consequence of the disappearance of the PN barrier l(22|) 
for small A; in the Ap^ — > limit, the restoring force 
vanishes for fluctuations in the direction parallel to the 
minimal energy path. 

In the d ~ A -1 <C 1 limit, the adiabatic correlations in 
the lattice dynamics at different sites disappear since the 
electron and the lattice deformation localize to a single 
lattice site. In the absence of adiabatic correlations in- 
volving multiple lattice sites there is no frequency soften- 
ing. Accordingly, in Fig. [3] one observes that for 1/A — > 
all the frequencies up approach the bare value luq. 



E. Curvature of the minimal energy path 



with 



«o(s) = - 



dQg(s) da(s) 
ds ds 



Us) 



d 2 d(s) 

ds 2 



(31) 



It is assumed here that the second term inside the 
square brackets of Eq. (|30j) is small for lattice configu- 
rations Q a relevant for the polaron dynamics. Other- 
wise, the singularities appear in the theory, indicating 
that the regime, for which the moving set of coordinates 
{s,Q a (s)} is inadequate, has been reached. In particu- 
lar, in the case when the fluctuations orthogonal to the 
minimal energy path are described by the zero-point mo- 
tion, (Q a ) ~ 1) ° ne finds that the applicability of the 
current formalism is restricted by the condition 



For the lattice with cyclic boundary conditions, the 
minimal energy path forms a closed loop in the configu- 
ration space of lattice deformations u. The curvature of 
this path, which depends nontrivially on e p and A, intro- 
duces a kinematic coupling^ between the motion along 
the minimal energy path and the motion orthogonal to 
it. This effect can be investigated by expressing the ki- 
netic part of the Hamiltonian lfT4|) in the representation 
of the moving set of coordinates {s, Q a (s)}. 

First, in the moving set of coordinates the lattice de- 
formation u is expressed as 



U(S, Q a ) = a(s) + 22 Qa(s) Qa 



(27) 



with a(s) the lattice deformation corresponding to the 

point s along the minimal energy path and Q a (s) the 
unit vector representing the unit displacement along the 
normal coordinate Q a , orthogonal to the minimal energy 
path. The total differential du is given by 



du = 



da(s) 
ds 



E 



dQg(s) 

ds 



Qa 



ds + Q a (s) dQ a 



(28) 

Multiplying Eq. (|28| by the unit vector da(s)/ds, which 
is tangential to the minimal energy path, gives 



da(s) 
ds 



du = 



1 



E 



dQa(s) da(s) 
ds ds 



Qa 



ds . (29) 



Using the completeness relation for the unit vectors 

{da(s) / ds,Qa\ , after some algebra, one obtains the fol- 
lowing expression for the kinetic energy associated with 
the motion along the minimal energy path§£ 



ds 2 



1 - E Ka ( s ) ( 3c 



(30) 



k(s)«1, k(s) = ,/E k «( s ) 



(32) 



where the coefficients K a (s) determine how the curvature 
of the minimal energy path is partitioned among the nor- 
mal modes a. 

Values of the coefficients k q (s) for different modes may 
be estimated by performing a simple analysis of the adi- 
abatic potential Uad(u). In this respect, it is convenient 
to use quantities which are simple functions of Hamilto- 
nian parameters. In particular, it is convenient to con- 
sider the relationship between two lattice deformations 
associated to the neighboring stationary points of the adi- 
abatic potential. The displacement h PN from the min- 
imum to the nearest saddle point in the configuration 




2 I 



Figure 4: (Color online) In general, the projection of the dis- 



placement vector h , given by Eq. l133[l . is substantial only 
in the direction of the pinning /3 = P and breather /3 — B 
coordinates (dashed curves), cos 2 jp + cos 2 jb ~ 1. For large 
polarons (A < 1) it is substantial only in the direction of the 
pinning coordinate, cos 7p ~ 1. The PN barrier for e p = 10 
is plotted by solid curves. 
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space u in terms of the normal modes, obtained in Eq. 
(|26)l . is given by 



PN 



tPN 



:pn< 



^2Qp COS 7/3 



(33) 



with \h PN \ ~ g/coo and the arguments 7/3 are functions 
of A. Since K a (s) ~ ug/g, it is the arguments 7/3 that 
define the nontrivial relation between the curvature of 
the minimal energy path and the normal modes f3. 

From Fig. [4] one can see that h PN lies almost entirely 
in the plane of the lattice configurations spanned by the 
pinning (3 = P and breather (3 = B coordinates. That 
is, cos 2 7/3=p + cos 2 7/3=5 > 0.99 is satisfied for all A. 
According to these findings, one may conclude that the 
effects due to the curvature of the minimal energy path 
are, in general, important just for the two normal modes 
P with the lowest frequencies. The interpretation of this 
result in the context of the moving set of coordinates is 
straightforward. For the kinematic coupling in Eq. (J30j) , 
only the coefficient n a (s) associated to the moving nor- 
mal mode a with the lowest frequency is significant, while 
for the other modes the coefficients can be neglected. 



V. LARGE ADIABATIC POLARONS 

The large adiabatic polaron (d/a ~ A -1 ^> 1) has 
many analogies with a classical soliton (with a funda- 
mental difference that it carries the elementary charge 
e). Its local dynamics is characterized by a large gap in 
the electron spectrum A ad ^> ujo, involving a large lat- 
tice deformation with respect to the zero point displace- 
ments. Concerning the translational dynamics, the quan- 
tum tunneling effects can be neglected because the pin- 
ning potential (|22f vanishes exponentially with d. Thus, 
for A ad 3> Wo and A -C 1, the same arguments as those 
appearing in the soliton theory, based on the large effec- 
tive mass, can be used to justify the use of the classical 
approximation. 

It should be stressed that long-range adiabatic corre- 
lations occur only for one dimensional Holstein polarons. 
Indeed, the regime of large adiabatic polarons makes the 
ID case exceptional. In contrast, the other regimes in 
the phase diagram discussed herein appear regardless of 
the dimensionality. 

In what follows, a short overview of the properties of 
the large adiabatic polaron derived in the classical limit 
is given in Sec. IV Al The crossover from large adiabatic 
polarons to large polarons in the weak-coupling regime 
is discussed in the context of the continuum (A <C 1) 
approximation in Sec. IV Bl The crossover from freely 
moving to pinned adiabatic (A ad 3> u>q) polarons due 
to the lattice coarsening effects is discussed in Sec. IV CI 

For both of the afore mentioned crossovers involving 
large polarons one should consider quantum fluctuations 



of the lattice deformation field, the nature of which is 
qualitatively different for the two cases. By crossing 
towards weak-coupling (Apjy ~ wo), it turns out that 
the nonadiabatic quantum electron-phonon correlations 
take control of the polaron dynamics. On the other 
hand, the polaron dynamics remains adiabatic in the 
crossover regime between freely moving and pinned po- 
larons (A <~ 1). In this latter case, the quantization of the 
lattice deformation field is necessary in order to describe 
the quantum tunneling through the PN barrier. 



A. Classical lattice theory of large adiabatic 
polarons 

By treating the lattice classically in the continuum 
limit for the polaron solution with the minimal energy, 
one obtain a 53 i 72 



= — \Vn\ 2 , Vn = ^ scch [A (n - £/a)/2] . (34) 

is the classical lattice deformation at the site n and 
T] n is the electron wave function. A defines the electron- 
phonon correlation length d, whereas £ represents the 
position of the polaron on the continuous axis. The 
binding energy of the large adiabatic polaron lj34j) is 
one third of the gap in the adiabatic electron spectrum, 
A po i = Aad/3, with A A d = £ P A/4IS. 

The solution in Eq. ([341 , which breaks the (full) trans- 
lational symmetry, defines the static polaron properties. 
The electron r\ n is trapped by the self-consistent lat- 
tice deformation u^ 1 . The full translational symmetry 
is restored by considering £ as a dynamical coordinate, 
£ = x — vt, where £ corresponds to the Goldstone mode 
of the continuous translational symmetry. In the contin- 
uum approximation, the effective polaron mass m po ^, 



m po i ~ (du n /d£) 
m po i 



m et 



2 1 ( u 2 



(35) 



The same 



is linear in the mass of nucleus, m po i ~ u 
result may be derived in the context of the moving set 
of coordinates {s,Q a (s)} by neglecting the curvature of 
the minimal energy path, 



ds = s [ds/d£] d£ 



ds 1 



(36) 



Furthermore, such behavior of the effective mass together 
with the behavior of the polaron size d ~ 1 / A is obtained 
from the scaling properties of the infinite-order diagram- 
matic perturbation theory^, which is translationally in- 
variant and quantum from the outset. 

Concerning the local dynamics, the elastic constants 
u>/3 and the parity of the normal modes derived in the 
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continuum approximation (Apjv = 0) are given in Table 
EP^. The pinning mode, for which the elastic constant 
uif3=p is zero, corresponds to the Goldstone mode for the 
polaron translation. The values in Table U represent the 
A" 1 -> limits of the plots in Fig. H 



B. Crossover towards weak-coupling regime 

The properties of the large adiabatic polaron have been 
discussed so far by treating the translational and local 
dynamics separately. However, due to the curvature of 
the minimum energy path, a kinematic coupling emerges 
between the coordinates s and Q a (s). Within the contin- 
uum approximation, the curvature of the minimal energy 
path dM]) scales as^ 



(37) 



By inserting Eq. (j37|) into Eq. (|32|) , the regime of weak 
kinematic coupling can be identified, 



< 1 



£p A > UJ ■ 



(38) 



One observes from Eq. f38|) that, in the continuum 
A <C 1 limit, the condition lj32]) is just a restatement of 
the condition Aad = £ p A/4 > wo that justifies the adi- 
abatic approximation. In other words, for large polarons 
the singularities in the kinetic energy ([30)) obtained in the 
adiabatic approximation coincide with the breakdown of 
the adiabatic approximation. This result has a simple 
physical interpretation. The quantum effects associated 
with the zero-point motion perpendicular to the minimal 
energy path should be a small perturbation to the adia- 
batic motion of the equilibrium lattice profile (|34| along 
the continuous axis. If this condition is not satisfied, 
the adiabatic approximation breaks (Aad ~ ^o) and the 
nonadiabatic quantum fluctuations take over. 

This interpretation is consistent with Eq. J9j, derived 
from the weak-coupling side. It can be therefore con- 
cluded that the crossover between large adiabatic po- 
larons and large nonadiabatic polarons in the weak- 
coupling regime occurs for e p A ~ ujq, or equivalently, 
for g/uo ~ (t/uj )i. 



Table I: The space symmetry and the frequencies of the adi- 
abatically softened normal modes obtained in the continuum 
approximation. 







UJ2 


U>3 




symmetry — 


+ 




+ 


+ 


U>/3 /tOO 


0.648 


0.882 


0.912 


0.949 0.958 



C. Crossover towards pinned polarons 

While in the continuum limit the large adiabatic po- 
laron behaves essentially as a free particle, the polaron 
spectrum exhibits a gapped band structure when the lat- 
tice coarsening effects are taken into account. In the sim- 
plest terms, the pinning effects on the polaron dispersion 
properties can be analyzed by neglecting the curvature 
of the minimal energy path for the polaron translation, 
which is appropriate in the regime of large adiabatic po- 
larons (e p » £ p A » uj ). 

Assuming that the PN potential behaves approxi- 
mately as Upn(s) = Apjy sin 2 (ns/s a ), with s a being 
the length of the minimal energy path corresponding to 
the translation of the polaron by one lattice site, 



m,p i 



(39) 



the polaron dispersion properties are determined by 
Mathieu's equation. Introducing x = s/s a , one obtains 



4 s 



1 d 2 

2 Q x 2 



E 



(n) 
K 



Apjy sin 2 (tt x) 



$>(x)=0. (40) 



The translationally invariant solutions of the Math- 
ieu's equation (|40|) . shown in Fig. [5j describe the disper- 
sion properties of large adiabatic polarons in the presence 
of pinning effects. The band structure in Fig. is ob- 
tained by calculating numerically the stationary points 
of Uad(u) that define s a in Eq. lf39j) and Apjy in Eq. 
(|2i"]) . In Eq. {20]), is the energy of the polaron with 
momentum K, with n distinguishing between different 
bands. 

In Fig. [U the PN barrier is plotted by the thick curve. 
For large PN barriers [right side of Fig. [5] the polaron 
bands are very narrow. In this limit, the separation be- 
tween different bands is given by the frequency of the 
pinning mode uj a= p, discussed in connection with Eq. 
(f26j . In the limit of a vanishing PN barrier [the left side 
of the Fig. [5] all the gaps between bands close. The 
spectrum describes a free-like particle, with the effective 
mass given by Eq. |39|) . In the continuum A <C 1 limit, 
Eq. (|39|l reduces to Eq. ([35|) . the latter being derived in 
the classical approximation. That is, in the continuum 
A « 1 limit, the same effective mass characterizes the 
quantum and the classical motion of the large adiabatic 
polaron. 

The band structure in Fig. [5] involves the bands for 
which all the oscillators associated with the dynamics 
orthogonal to the energy path are in the ground state. 
Additional bands are obtained by considering the excita- 
tions of these oscillators, 



E(K,n,n a ) = 



a^P 



(41) 
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Figure 5: (Color online) The translationally invariant solu- 
tions of Eq. iHOl) that describe the dispersion E)£' of large 
adiabatic polarons in the presence of the PN barrier. K is 
the momentum and n is the band number. For n even, where 
n = denotes the lowest band, the bands are bounded from 
below and above by the K = (full curves) and K = tt 
(dashed curves) states, respectively. For n odd the situation 
is opposite, i.e., the bands are inverted. The thick curve is 
the PN barrier Apat. 

The treatment of the large adiabatic polaron in Eq. lj4lj) 
neglects entirely the effects due to the curvature k(s) 
of the minimal energy path. With these effects taken 
into account, one would observe a weak hybridization 
between different n a ^ bands. The hybridization issue 
is discussed in more detail in Sec. lYIBl 



VI. SMALL PINNED POLARONS 

A. Adiabatic translational dynamics 

It was argued in Sec. IV Bl that, in the context of the 
(adiabatic) moving set of coordinates, the increasing im- 
portance of kinematic contributions in the continuum 
limit can be traced from the curvature of the minimal en- 
ergy path k(s). In particular, the singularities in Eq. f30|) 
are avoided only for sufficiently large e p , s p /ojo > A -1 . 
For small polarons, on the other hand, the singularities 
in the context of the moving set of coordinates appear 
due to the particular shape of the adiabatic potential for 
any value of e p . That is, although in the e p 3> limit 
k(s) satisfies Eq. |32|) along the minimal energy path 
that connects two minima of the adiabatic potential, this 
quantity is singular at the minima. As argued here, the 
nature of the singularities in the cases of the large and 
small polarons is fundamentally different. In the former 
case, the quantum fluctuations perpendicular to the min- 
imal energy path are involved, whereas in the latter the 
singularities are found along the minimal energy path it- 
self. 

The curvature k(s) is singular at the minima of Uad(u) 




A,= 1.2 




X=4 



Figure 6: The adiabatic potential as function of coordinates 
Q/3=p and Qf3=B, belonging to the pinning and breather nor- 
mal modes for A = 1.2 and A = 4. The absolute mini- 
mum of Uad{Qi3=p,Qi3=b) is at the origin, while the two 
local minima correspond to polarons centered at the nearest- 
neighboring sites. 



when the normal mode with the lowest frequency is the 
(even) breather mode. As can be observed from Fig. 
[21 this happens for A > 1.4. The shape of the adiabatic 
potential as a function of the pinning Qp=p and breather 
Q/3=b normal coordinates is compared for large {uip=p < 
W/3=b) and small (utp—p > lop=b) polarons in Fig. [6l In 
Fig. [6^ the minimal energy path passes the minimum 
u M tangentially in the direction of the axis Q/3=p. On 
the other hand, in Fig. [6Jd the minimal energy path in 
the vicinity of u M is in the direction of the axis Qp = B, 
with u M representing the (singular) turning point. This 
result regarding the appearance of singularities in the 
small polaron case can easily be generalized to higher 
dimensional lattices. 



B. Crossover between small and large adiabatic 
polarons 

With the complex shape of the adiabatic potential, 
the quantitatively accurate description of the adiabatic 
translational dynamics for small polarons relies on nu- 
merical methods. The problem is simplified by the 
large PN barrier, which permits one to consider only the 
nearest-neighbor hopping processes. One may start with 
a set of vibrational wave functions that describe har- 
monic fluctuations around the minima of Uad(u), and 
then calculate the overlap integrals between the vibra- 
tional wave functions involving the nearest neighbor sites. 
Because the minimal energy path within the unit cell 
practically lies in the plane of the pinning and breather 
modes (as discussed in connection with Fig. [4]), the over- 
lap integrals should be considered explicitly only for these 
two modest. Following such a procedure, one obtains a 
tight-binding problem that can easily be solved, in partic- 
ular numerically. In fact, with generalizations necessary 
for the description of large polarons and nonadiabatic 
contributions, the correspondingly extended approach is 
used by the relevant coherent states method (RCSM). 
The latter methodology, introduced recently in Refi^, 
calculates accurately the polaron spectra for the whole 
parameter space. 

Figure [7K shows the RCSM polaron spectrum in the 
crossover regime Awl between large and small adia- 
batic polarons (t/u>o — 250). It involves the lowest po- 
laron band, denoted by G, and six excited bands. All 
the states are given with respect to the polaron ground 
state Eq = E^T^ . The notations B (breather) and P 
(pinning) relate the narrow bands on the right side of 
the spectra (large PN barrier) with the harmonic adia- 
batic theory. For example, PB denotes the band which 
involves the excitation of the pinning and the breather 
mode shown in Fig. [3J i.e., its position in the spectrum 
is determined by u>p = B + wp = p. BB means a double 
excitation of the breather mode, and so on. 

The spectrum in Fig. [7K is given for the same choice of 
parameters as in Fig. [H By comparing these two spec- 
tra for A < 1, one finds that the energies of the (large) 
polaron states in Fig. [7^ behave according to the pre- 
diction of the approximate calculation of Eq. l(4T|l . de- 
rived by entirely neglecting the curvature of the minimal 
energy path. In particular, there is no significant dif- 
ference in the dispersion properties of the large polaron 
described by the solutions [Fig. [5] of Mathieu's equa- 
tion (|40|l and the RCSM states marked by the arrows 
in Fig. [7^. Concerning the local dynamics, the K = 
state marked by the filled circle involves the excitation 
of the breather mode. Its energy on the left side of Fig. 
EH (A < 1) is very close to the energy of the breather 
mode iop=B ~ 0.65, obtained in the continuum approxi- 
mation [see Table [J. This result agrees with Eq. l(4lj) as 
well. Thus, it can be concluded that our Eq. ([41]) gives 
a very satisfactory description of the large polaron band 
structure in the presence of lattice coarsening effects. 
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Figure 7: (Color online) The lowest polaron band and the 
excited polaron bands associated to the pinning and breather 
modes, calculated by the RCSM, for t = 250 in the upper and 
for t = 25 in the lower panel (ujo = 1). The spectra are shifted 
by the ground state energy. The K = and K — tv states 
are plotted by full and dashed curves, respectively. The inset 
shows the band structure for t = 25 in terms of 22 states with 
different momenta for each band, K = 0,ir,n x 0.15, n < 20. 



The polaron spectrum in the crossover regime between 
pinned small polaron states and freely moving large po- 
laron states exhibits a particularly complex structure. It 
is characterized by a strong hybridization between dif- 
ferent excited bands, as can be seen from Fig. [7b for 
t/ujo = 25. The hybridization is best observed from the 
inset of Fig. [7b, in which the each band is shown in terms 
of 22 states with different momenta. The first and second 
excited bands touch in Fig. [7b for g/too ~ 6.3. This band 
touching corresponds to the symmetry allowed crossing 
of the two K = it polaron states of opposite parity. A 
detailed discussion of parity properties can be found in 
Ref.— . For g/ujQ < 6.3 in Fig. [7b, the strong hybridiza- 
tion results in a very atypical dispersion; the K = and 
K = 7r states of the first excited band intersect twice. 

On the left side of the spectrum in Fig. [7b involving 
large polarons (t/u>o = 25), the excited band correspond- 
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ing to the excited breather mode (the K = state is 
denoted by the filled circle) shifts towards the phonon 
threshold with decreasing g. This indicates the presence 
of nonadiabatic contributions to the polaron dynamics, 
the consequence of which is a suppression of the adiabatic 
softening of the breather mode that moves with the po- 
laron. As one can see from Fig. [7K, this shift is far less 
pronounced for larger values of Aad/^o [t/u>o = 250 in 
Fig. vs. t/ujQ = 25 in Fig. EH- 

C. Adiabatic and nonadiabatic small-polaron 
hopping 

In the regime of small pinned polarons e p A > e p » ui , 
the large gap in the adiabatic electron spectrum, Aad ~ 
e p , makes the local dynamics adiabatic irrespectively of 
the ratio t/uj . However, due to the large PN barrier, 
Apn ~ Aad/4, two kind of processes should be consid- 
ered for the translational dynamics. The small polaron 
hopping is adiabatic when dominated by the processes in 
which the electron and the lattice deformation move to- 
gether to the neighboring site. Yet, when the deformation 
field is pinned by the PN barrier too strongly, the dereal- 
ization energy is gained by the nonadiabatic rather than 
the adiabatic hops of the electron to the neighboring site. 
During the nonadiabatic hops, the electron is detached 
temporarily from the lattice deformation field. 

It is convenient to discuss the distinction between the 
adiabatic and nonadiabatic regimes for the small-polaron 
hopping by considering the two-site model as the results 
exist in closed for m 75 ! 76 ' 77 ' 78 . As long as the electron- 
phonon correlation length is shorter than the size of the 
cluster under consideration, one may expect that the 
obtained solutions satisfactorily reproduce the main po- 
laron properties. In particular, in the limit d ~ 1/A — > 0, 
it can be argued that the dimensionality and the cyclic 
boundary conditions that characterize clusters of differ- 
ent size are not of decisive importance. 

Regarding the nonadiabatic small-polaron hopping, it 
should be noticed that the hopping energy in the lead- 
ing order in t, given by Eq. (fit)]) , depends only on the 
quasi-particle weight exp (—e p /uio). All the phonons con- 
tributing to this quantity involve a single lattice site, 
which makes it independent of the cluster size and the 
dimensionality of the lattice. 

By separating out the homogenous lattice deformation, 
the adiabatic potential for the two-site cluster is obtained 
as2£ 



U 2 (x) = e p (x 2 - v / 25 2 + A- 2 ) , x = uj x/g , (42) 

with x the relative lattice deformation between the two 
sites. Obviously, the effects due to the curvature of the 
minimal energy path cannot be addressed within the two- 
site model. However, for A 3> 1, the adiabatic hopping 
energy that characterizes the motion along the path that 
connects two minima behaves similarly, independently of 



the cluster size and system dimensionality. Namely, the 
scaling properties of the adiabatic potential, given by Eq. 
(fTF)) . are a general property of the Holstein Hamiltonian. 
For example, the adiabatic potential for the three-site 
cluster in polar coordinates takes the formZ2: 



U z {f,(t>) = 2e p y2V3cos(0/3) 

f = uj r/g , coscf) = -f/h% , 
h = f 2 +3A~ 2 /2 

f = f 3 cos(30) + (3A- 2 /2) i ■ 

For both, U 2 (x) and U^{f, <j>), only A defines the essential 
dependencies. Starting with the 1/A — > limit, in which 
the PN barrier is the same for any cluster, Apjy = s p /2, 
one finds that the leading correction in 1/ A involves the 
top of the PN barrier. Namely, at this point the polaron 
is centered between two lattice sites, meaning that the 
electron is equally shared between the two neighboring 
sites. The kinetic energy gain due to such electron der- 
ealization is given by the electron hopping energy t. On 
the other hand, it is easy to verify that the 1/A correc- 
tions to the minima of the adiabatic potential are given in 
the second order. Thus, the PN barrier up to the leading 
order in 1/ A behaves according toApjy ~e P (1/2 — 1/A), 
in perfect agreement with numerical results in Fig. ((2|), 
derived in the context of the infinite lattice. 

For A > 1 (pinned polarons) the potential U 2 (x) ex- 
hibits two minima at x M = ±\/I — \~ 2 /V2, which are 
separated by the PN barrier. The frequency of the har- 
monic vibrations around the two minima is softened with 
respect to the bare frequency ojq, 



uj'/ujo « y/l-X' 2 . (43) 

The polaron hopping energy can be estimated by using 
the quasi-classical WKB approximation^, 



ln(tf ol ) ~ -2 / V\Eo-U 2 (x)\ dx . (44) 
Jo 

with E = U 2 {x M ) + ui'/2. Following Ref3, for g > t 
the adiabatic hopping energy is obtained as 

t% « (27r)-Vexp(-( ep /wo)/(A)), 
/(A) = v / T T A^-A- 2 ln^ 1 + V A 1 _7 A '^ .(45) 

In the A 1 limit, which is of interest here, t pol scales 
as 

Hto' 

tpoi * — exp (-E p /hu (1 - 1/2A 2 )) (46) 
~ exp(-e p /uj ) cxp(t 2 /2.g 2 ) . (47) 
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The term t 2 /g 2 in the exponential does not appears in 
the case of nonadiabatic hopping in Eq. (fit)]) . Thus, it 
follows that for t > g the adiabatic hopping dominates 
the small-polaron translational dynamics&Si. 

On the other hand, for g > t, the integration in 
Eq. (|44|) gives the same exponential behavior as in 
Eq. 



t 



pol 



g exp (-Sp/wo) 



(48) 



In the regime g >t, when the adiabatic and nonadiabatic 
hopping energies differ only in the preexponential factor, 
the nature of the small-polaron hopping is apparently 
mixed. Yet, since Eq. (|48"1) is derived by neglecting kinetic 
energy contributions related to nonadiabatic effects, its 
applicability in the whole g > t regime still remains to 
be verified. 

As pointed out in connection with Eq. I|14p . the lattice 
kinetic energy operator in the adiabatic limit behaves as 
if it commutes with the electron wave function. In gen- 
eral, of course, this is not the case and one should con- 
sider the corresponding kinetic contributions as well. In 
particular, for the two-site model the calculations can be 
performed analytically. The ground-state electron wave 
function corresponding to the adiabatic potential ([42]) is 
given byZ£ 



|?/) = cos rj(x) \c\) + sin rj(x) |c|) 



cos rj(x) 
sin 77(f) 



V2x 



1 



2 2V2Z 2 + A" 2 



V2x 



1 



2 2\/2~PTA~ 



(49) 



where |c{) and \c\) are the wave functions of the electron 
localized at the site 1 and 2, respectively. The exact 
expectation value of the lattice kinetic energy is obtained 
as 



(V\T\V) 
AT(5) 



w 



A" 



A- 2 + 2x 



(50) 



The nonadiabatic contribution AT(i) has the shape of a 
squared Lorentzian centered at the top of the PN barrier 
x = 0. Using a simple variational approach, it is shown in 
Refi^ that AT(i) introduces corrections to the adiabatic 
hopping energy (|48| given by 



9 e 



TTOJQ 
8 t 



O 



(51) 



The uJo/t correction in Eq. f5"Tj) clearly demonstrates 
that for t < uj the adiabatic treatment of the small- 
polaron hopping breaks down. In this regime the function 
^(cc) + AT(i) at x w takes values that are larger 
than the free electron energy, meaning that the electron 
necessarily detaches from the lattice deformation during 
the hop. 

The nature of the small polaron hopping is controver- 
sial in the literature^ 2 -. From the current analysis, it 
might be concluded that the translational dynamics of 
the small pinned polarons (e p A > e p ^> w ) is nonadia- 
batic for t < u>q and adiabatic for t > g. The regime 
wo t < g appears as a crossover regime in which 
the hopping processes involve a mixture of adiabatic and 
nonadiabatic contributions. The situation is much sim- 
pler with the local dynamics. In the whole regime of small 
pinned polarons the nonadiabatic contribution AT(x) 
given by Eq. lj50j) vanishes for x « x , i.e., the local dy- 
namics is adiabatic. In particular, for the two-site model 
the local dynamical correlations are characterized by the 
adiabatically softened frequency (|43|) . 



D. Polaron self-trapping 

The term polaron self-trapping is frequently used in the 
literature to mark the change from free-particle-like to 
exponentially-large effective mass behavior of the polaron 
states. In the small t limit, t/ujo <C 1, discussed in Sec. 
IIIIB1 this change occurs for a nonadiabatic combination 
of parameters, g/u>o ~ 1. For the opposite limit, t/u)o^> 
1, it is argued in Sec. IVI Bl that the polaron self-trapping 
occurs for Awl, where A is a parameter independent of 
the mass of the nucleus. 

By constructing a very simple additive combination 
of the nonadiabatic and the adiabatic criteria for the po- 
laron self-trapping, one obtains an empirical formula that 
accurately describes the polaron self-trapping regime for 
the whole parameter space, 



gsr 



(52) 



The formula ([52)1 was suggested in Ref4^, by considera- 
tion of numerical solutions for the lowest polaron band. 
It satisfies both limiting behaviors discussed in this work: 
gsr/ wo ~ 1 in the small t limit and g% T /t loq = Xst ~ 1 
in the large t limit. 

For moderate values of t/uo, for which the two terms 
contributing to gsr in Eq. i|52|) are comparable, a com- 
plex mixture of adiabatic and nonadiabatic contributions 
takes place in the polaron self-trapping regime. In Fig. 
[5] the polaron spectrum is shown for t/tJo = 2.5. Unlike 
in Figs. [IK and [TJd, the gap between the lowest and the 
excited bands remains considerable for all couplings, in- 
dicating that the adiabatic dynamics is characterized by 
a substantial PN barrier. That is, whereas this gap closes 
in Fig. [7K as g is decreased [see the lowest two K = it 
states pointed by the arrow], in Fig. [8] all the excited 
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Figure 8: (Color online) The polaron band structure for t = 
2.5 (uio = 1) calculated by the RCSM. The spectrum is shifted 
by the ground state energy. The K = and K = it states are 
plotted by full and dashed curves, respectively. 

polaron bands shift towards the phonon threshold. This 
shift is a clear indication that the nonadiabatic correla- 
tions take control of the local dynamics in the left-part of 
Fig. meaning that this parameter regime corresponds 
to weak couplings. 

On the right side of Fig. [8l corresponding to the pinned 
small polarons, the positions of the narrow excited bands 
in the spectrum are well predicted by the harmonic adi- 
abatic theory. Since t ~ g, according to the discussion in 
Sec. IVI C\ the polaron hopping energy involves adiabatic 
and nonadiabatic processes. This conclusion agrees with 
the numerical analysis of the small polaron hopping for 
the infinite lattice carried out in Refi^. 



VII. PHASE DIAGRAM 
A. One dimensional case 

The discussion of the polaron properties in this work 
can be summarized in terms of the two dichotomies in the 
polaron dynamics, local vs. translational and adiabatic 
vs. nonadiabatic, with emphasis given to the importance 
of the lattice coarsening. The corresponding phase dia- 
gram for the Holstein polaron is proposed in Fig. [9l Dif- 
ferent regimes and crossovers between them are d rawn in 
the two-dimensional parameter space g/uo vs. y^t/wo- 

Along the solid curve in Fig. [9J corresponding to 
Eq. ([9}, the polaron binding energy A po i is of the order 
of the bare phonon energy wo- Above this curve the lo- 
cal dynamics is adiabatic (self-trapped electron) , whereas 
below is the weak-coupling regime (weakly dressed elec- 
tron). It is stressed in Sec. IIII Al that. for weak couplings, 
the local and translational dynamics of the polaron are 
nonadiabatic. In the absence of adiabatic correlations, 
the spectrum below the phonon threshold exhibits a sin- 



gle polaron band. Additional polaron bands below the 
phonon threshold appear in the regime above the solid 
curve in Fig. [9j due to the adiabatic softening of the 
phonon modes associated with the (large) lattice defor- 
mation. 

Along the thin dashed curves in Fig. [9] the electron- 
phonon correlation length d (polaron size) takes constant 
values. In the weak-coupling regime d is independent of g, 
whereas in the regime above the solid curve d is indepen- 
dent of the lattice mass, being constant for X = g 2 /ojq t 
constant. The thin dashed curves in Fig. [9] are given for 
A = 2™, with n representing an integer in the range from 
—2 to 7, inclusively. Along the thick dashed curve, corre- 
sponding to A = 1 , the polaron deformation spreads over 
a few lattice sites. The polarons are small/large in the 
region on the left/right side of the A = 1 curve. In the 
case of large polarons, as argued in Sec. IIII Al for weak 
couplings and in Sec. |V] for the adiabatic limit, in order 
to calculate the polaron binding energy and the effec- 
tive mass one may apply the continuum approximation, 
which neglects the lattice coarsening effects. 

The green area in Fig. [9] corresponds to the regime of 
large adiabatic polarons. In this regime, the polarons be- 
have as free particles with large mass. The mass is given 
by Eq. derived here by treating the lattice classi- 

cally in the continuum approximation. In the regime of 
large adiabatic polarons, as discussed in connection with 
Eq. ([38]) , the kinematic coupling in Eq. (|30|) between the 
translational degree of freedom and the adiabatically- 
softened normal mode s that move with the polaron is 
weak, scaling as \/A _1 wo/.g- This means that in the case 
of large adiabatic polarons, the local and translational 
dynamics can be treated separately, as in Eq. (|4Tjl . From 
Fig. [9j one observes that the large adiabatic polaron 
crosses directly into a large nonadiabatic polaron for 
weak-couplings (nonadiabatically dressed electron). The 
position of this crossover in the phase diagram follows 
from Eqs. ((9} and (|38| . i.e., it is given by t/u>o ~ {g/uio)* . 

According to the nature of the translational dynamics, 
the regime of small pinned polarons (e p X > e p ^> ujq) in 
Fig. [9] is divided by the dot-dashed and dotted curves, 
corresponding to t ~ wo and t ~ g, respectively. As ar- 
gued in Sec. IVI C| for t > g, the small-polaron hopping is 
dominated by adiabatic processes, in which the electron 
and the lattice deformation hop together to the neighbor- 
ing site. For t < ojq, the nonadiabatic hopping prevails 
because the electron is too slow to follow the lattice de- 
formation during the hop. For luq < t < g, the polaron 
hopping involves a mixture of adiabatic and nonadiabatic 
processes. 

The circles in Fig. [9] are drawn along the curve defined 
by Eq. ([52"ll . The circles mark the polaron self-trapping, 
i.e., the crossover regime between pinned polarons (above 
the circles) and freely moving polarons (below the cir- 
cles). In the adiabatic limit t 3> u>q, the polaron self- 
trapping corresponds to the large circles, drawn close to 
the A ~ 1 curve. For A < 1, the large adiabatic polaron is 
free to move since the PN barrier is ineffective, decreas- 
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Figure 9: (Color) The phase diagram for the ID Holstein polaron plotted in the parameter space g vs. yi, with ljo as the 
unit of energy. The weak-coupling regime (lattice quantum fluctuations important) [Sec. IIII Aj . the regime of large adiabatic 
polarons (classical behavior) [Sec. [V], the regime of small adiabatic polarons (adiabatic hopping) [Sec. I VI Aj . the regime 
of small nonadiabatic Lang-Firsov polarons (nonadiabatic hopping) [Sec. IIII Bj are shown in blue, green, yellow and red, 
respectively. The regime of small polarons [Sec. IVICj . for which neither the adiabatic nor nonadiabatic contributions prevail in 
the translational dynamics, corresponds to the transition area from red to yellow. Along the thick and thin dashed curves the 
electron-phonon correlation length d (polaron size) takes constant values. For weak couplings d is independent of g, whereas 
for strong couplings (adiabatic local dynamics) d is constant for A = g 2 /ujq t constant. The thick dashed curve denotes the 
crossover between large (A < 1) polarons (continuum approximation applicable) and small (A > 1) polarons. The nonadiabatic 
to adiabatic crossover in the local dynamics is represented by the solid curve, g 2 ~ 1 + 2\/t. This crossover is derived from the 
weak-coupling side in Eq. ([9]) and from the strong coupling side in Eqs. I|38l) and (lllll for large and small polarons, respectively. 
The circles denote the polaron self-trapping g ~ 1 + \ft [Eq. (I52II ]. marking the crossover between pinned polarons (above 
the circles) and freely moving polarons (below the circles). The two crossovers involving small pinned polarons denoted by 
the dotted (g ~ t) and dot-dashed (t ~ 1) curves are associated with the changes in the translational dynamics. The small 
polaron hopping is nonadiabatic for t < 1 and adiabatic for t > g, while for 1 < t < g the hopping processes involve a mixture 
of adiabatic and nonadiabatic contributions [Sec. IVI Cj . The white area corresponds to the regime of parameters in which 
polaron formation involves a complex mixture of adiabatic and nonadiabatic contributions and lattice coarsening effects. 



ing exponentially with the polaron size d ~ 1/ A [see Eq. 
(f22|l ] . For A <~ 1, as shown by Fig. [U the gapped polaron 
band structure develops from the free-particle-like spec- 
trum due to the pinning effects. On the other hand, for 
small polarons (A > 1), a strong hybridization between 
the excited bands associated with the lowest odd (pin- 
ning) and even (breather) modes occurs. As explained 
in connection with Figs. [4]and[6l the origin of this hy- 
bridization is the lattice coarsening. Namely, the coars- 
ening introduces not only the PN barrier in the polaron 
adiabatic motion along the lattice, but it also effectively 
couples the lowest two normal modes of the adiabatic 
lattice deformation (the pinning and the breather mode) 



through the translational motion. 

For t < loq, the polaron self-trapping is marked in 
Fig. [9] by small circles, corresponding to g rts (Jq. This 
crossover is identified in Sec. IIII Bl as the crossover be- 
tween the weak-coupling regime and the regime of small 
nonadiabatic (Lang-Firsov) polarons. More precisely, for 
t < ujo, the lattice pinning of the polaron coincides with 
the adiabatic self-trapping of the electron. The local cor- 
relations are adiabatic and nonadiabatic on the strong- 
coupling and weak-coupling sides of this crossover, re- 
spectively. On the other hand, on both sides of the 
crossover, the translational dynamics is characterized by 
the nonadiabatic polaron hopping energy, given by Eq. 



17 



crj. 

It is argued in Sec. I VI PI that the subtle balance be- 
tween nonadiabatic and adiabatic contributions on the 
discrete lattice describes the polaron self-trapping for 
moderate values of t/u/Q, i.e, 1 < t/u/Q < 20. In this 
regime of parameters, corresponding in Fig. [9] to the 
white area, the polarons exhibit particularly interest- 
ing properties. That is, although the static proper- 
ties, like the binding energy, change slowly with pa- 
rameters, the dispersion changes dramatically. For ex- 
ample, for t/u>o = 5, by increasing the coupling from 
g/uo = 3 to g/ujQ = 3.5, the polaron effective mass 
changes by two orders of magnitude, m po i/m e i f=a 0.147 
vs. m po i/m e i « 0.0019. At the same time, the change 
in the binding energy is far less dramatic, A po i/ujq ~ 2.5 
vs. and A po ;/wo ~ 4.5. 

B. Generalization to higher dimensional cases 

In contrast to the two and higher dimensional real- 
izations of the Holstein polaron, the particularity of the 
one dimensional D = 1 phase diagram is that it exhibits 
the regime of large adiabatic polarons. Namely, in the 
case of short range interactions, for dimensions greater 
than one the adiabatic polarons are unstable when the 
size of the polaron exceeds a critical value that depends 
on - D 69 ] 84 ] 85 , In particular for the Holstein model, for 
D = 2 one finds that the adiabatic polaron is unstable 
for A < 3.34 and for D = 3 it is unstable for A < 5.42^. 
In other words, for D > 1 the adiabatic polaron is stable 
only if it is small, i.e., pinned by the discrete lattice. 

On the other hand, for weak couplings, the polarons 
are large for t/wo 3> 1, irrespectively of the dimension of 
the system. By considering the generalization of Eq. ([9]) 
to arbitrary D^, 

e p /u < (t/oj ) D/2 , (53) 

one obtains, as a function of D, the part of the phase dia- 
gram in Fig. [9]that belongs to the weak-coupling regime. 
One observes that for D > 1 the large polaron for weak- 
couplings crosses into a small adiabatic polaron. 

Regarding small polarons, it is plausible to think that 
the structure of the corresponding part of the phase di- 
agram is unaffected by the dimension of the system. 
Namely, in the case of small polarons, the dispersion 
properties are determined by the nearest neighbor po- 
laron hopping energy, the nature of which cannot be af- 
fected fundamentally by the lattice dimensionality. In 
particular, both, the nonadiabatic and the adiabatic 
regime of small pinned polarons are reported in the con- 



text of the dynamical mean field theor y 86 ' 8 , which is 
exact in the infinite dimensional limit. 

VIII. SUMMARY 

The properties of the Holstein polaron are discussed in 
terms of limiting analytical results, while the recently in- 
troduced relevant coherent state method is used to gain 
detailed insights into the nontrivial mixing of fundamen- 
tally different contributions to the polaron dynamics in 
the crossover regimes (quantum vs. classical, weak vs. 
strong coupling, adiabatic vs. nonadiabatic, itinerant vs. 
self-trapped polarons, large vs. small polarons). New 
results are derived, particularly concerning adiabatic as- 
pects of the polaron properties. An original and unifying 
interpretation of the polaron formation is proposed with 
particular emphasis on lattice coarsening effects in terms 
of two dichotomies in the polaron dynamics, local vs. 
translational and adiabatic vs. nonadiabatic. The de- 
gree of applicability of various approximate approaches 
to the regimes encountered in the phase diagram is fully 
clarified, which gives a complete explanation of the low- 
frequency polaron band structure. 

For the one dimensional D = 1 Holstein model five 
regimes are identified and positioned in the parameter 
space. In the weak-coupling regime the local and the 
translational polaron dynamics are nonadiabatic due to 
quantum fluctuations of the lattice deformation field. In 
the regime of large adiabatic polarons, the Holstein po- 
laron has many analogies with classical solitons, moving 
freely along the lattice with a large effective mass. As 
a consequence of the discreteness of the lattice defor- 
mation field in the regime of small adiabatic polarons, 
the translational dynamics involves (adiabatic) quantum 
tunneling through the Peierls-Nabbaro (PN) barrier. In 
the regime of small nonadiabatic (Lang-Firsov) polarons, 
the translational dynamics involves nonadiabatic polaron 
hopping and adiabatic local correlations. Between small 
adiabatic and Lang-Firsov polarons, the regime of small 
pinned polarons is located, for which neither the adia- 
batic nor nonadiabatic contributions prevail in the trans- 
lational dynamics. For higher dimensional D > 1 cases, 
the large polarons exists only for weak couplings. As the 
coupling is increased, such large polarons cross directly 
into small adiabatic polarons. 
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